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We present an overview of the correlation-matrix methods developed recently by the CSSM Lat- 
tice Collaboration for the isolation of excited states of the nucleon. Of particular interest is the 
first positive-parity excited-state of the nucleon known as the Roper resonance. Using eigenvec- 
tors of the correlation matrix we construct parity and eigenstate projected correlation functions 
which are analysed using standardized methods. The robust nature of this approach for extracting 
the eigenstate energies is presented. We report the importance of using a variety of source and 
sink smearings in achieving this. Ultimately the independence of the eigenstate energies from the 
interpolator basis is demonstrated. In particular we consider 4x4 correlation matrices built from 
a variety of interpolators and smearing levels. Using FLIC fermions to access the light quark 
mass regime, we explore the curvature encountered in the energy of the states as the chiral limit 
is approached. We report a low-lying Roper state contrasting earlier results using correlation ma- 
trices. To the best of our knowledge, this is the first time a low-lying Roper resonance has been 
found using correlation matrix methods. Finally, we present our results in the context of the Roper 
results reported by other groups. 
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1. Introduction 

The first positive-parity excited state of the nucleon, known as the Roper resonance, (1440 
MeV) Pii, has been a long-standing puzzle since its discovery in the 1960's due to its lower mass 
compared to the adjacent negative parity, N"^ (1535 MeV) Sn, state. In constituent quark models 
with harmonic oscillator potentials, the lowest-lying odd parity state naturally occurs below the Pn 
state (with principal quantum number N = 2) ||l], 0], whereas, in nature the Roper resonance is 
almost 100 MeV below the Sn state. 

Lattice QCD is very successful in computing many properties of hadrons from first principles. 
In particular, in hadron spectroscopy, the ground states of the hadron spectrum are now well un- 
derstood. However, the excited states still prove a significant challenge. The first detailed analysis 
of the positive parity excitation of the nucleon was performed in Ref. [^] using Wilson fermions 
and an operator product expansion spectral ansatz. Since then several attempts have been made to 
address these issues in the lattice framework, but in many cases no potential identification of the 
Roper state has been made. Recently, however, in the analysis of Refs. [Q, ^, ^ a low-lying Roper 
state has been identified using Bayesian techniques. 

Another state-of-the-art approach in hadron spectroscopy is the 'variational method' ||7|, |8|], 
which is based on a correlation matrix analysis. The identification of the Roper state with this 
method wasn't successful in the past. However, very recently, in Ref. a low-lying Roper state 
has been identified with this approach employing a diverse range of smeared-smeared correlation 
functions. 

Here we discuss the new correlation matrix construction for isolating the puzzling Roper state 
and present our Roper results in the context of the previous results reported by other groups in 
recent times. 



2. Variational Method 

The two point correlation function matrix for p = can be written as, 

Gij{t) = (£Tr,p{r±(Q|;^,-WZi(0)|n)}), (2.1) 

X 

= Y,XfXf (2.2) 

a 

where, Dirac indices are implicit. Here, A" and A" are the couplings of interpolators Xi Xj 
at the sink and source respectively and a enumerates the energy eigenstates with mass m^. r± 
projects the parity of the eigenstates. 

Since the only t dependence comes from the exponential term, one can seek a Unear superposition 
of interpolators, ^yWy , such that, 

Gij{t + AO M« = e-'""^' Gij{t) uj, (2.3) 



for sufficiently large t and t + At. More detail can be found in Refs. [ [TOl , |rT| , [T2| ]. Multiplying the 
above equation by [Gjj{t)]^^ from the left leads to an eigenvalue equation, 

[{G{t))-'G{t + At)]iju'} = c''uf, (2.4) 
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where c"" = e is the eigenvalue. Similar to Eq.(2.4), one can also solve the left eigenvalue 
equation to recover the v" eigenvector, 

vf[G{t + At){G{t))-']ij = (2.5) 

The vectors and vf diagonalize the correlation matrix at time t and f + Af making the projected 
correlation matrix, 



vfG,y(0"?°<5«'^. (2.6) 



The parity projected, eigenstate projected correlator, 

Gl^vfGfj{t)u'j, (2.7) 
is then analyzed using standard techniques to obtain masses of different states. 

3. Lattice Details 

Our lattice ensemble consists of 200 quenched configurations with a lattice volume of 16^ x 



32. Gauge field configurations are generated using the DBW2 gauge action [ |13| , |14| ] and an 
(a) -improved FLIC fermion action [ |T5| ] is used to generate quark propagators. The lattice 
spacing is a = 0.1273 fm, as determined by the static quark potential, with the scale set us- 



ing the Sommer scale, ro = 0.49 fm [16|. In the irrelevant operators of the fermion action 



we apply four sweeps of stout-link smearing. Various sweeps (1, 3, 7, 12, 16, 26, 35, 48) of 



gauge invariant Gaussian smearing [17] are applied at the source (at t = A) and at the sink. 
The analysis is performed on ten different quark masses corresponding to pion masses mn = 
{0.797, 0.729, 0.641 , 0.541 , 0.430, 0.380, 0.327, 0.295, 0.249, 0.224} GeV. Error analysis is performed 
using a second-order single elimination jackknife method, where the ;^^/dof is obtained via a co- 



variance matrix analysis method. Our fitting method is discussed extensively in Ref. [ ]11[ ] 
The nucleon interpolator we consider is the local scalar-diquark interpolator [^,18], 



Xi{x) = e"^''{u'"{x)Cysd\x))u'{x). (3.1) 



4. Results 



We consider several 4x4 correlation matrices. Each matrix is constructed with different sets 
of correlation functions, each set element corresponding to a different number of sweeps of gauge 
invariant Gaussian smearing at the source and sink of the XiXi correlators. This provides a large 
basis of operators with a wide range of overlaps among energy states. 

We consider seven smearing combinations { 1=(1,7, 16,35), 2=(3,7, 16,35), 3=(1, 12,26,48), 



4=(3, 12,26,35), 5=(3,12,26,48), 6=(12,16,26,35), 7=(7, 16,35,48)} of 4 x 4 matrices. In Ref . [y] 
it was shown that one cannot isolate a low-lying excited eigenstate using a single fixed-size source 
smearing. The superposition of states manifested itself as a smearing dependence of the effective 
mass. 

In Fig.|I|, masses extracted from all the combinations of 4 x 4 matrices (from 1**' to 7"^) are 
shown for the pion mass of 797 MeV. Some dependence of the excited states on smearing sweep 
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Figure 1: (Color online). Masses of the nucleon, A^2 - states, from projected correlation functions as shown 



in Eq.(2.7) for the pion mass of 797 MeV. Numbers in the horizontal scale correspond to each combination 
of smeared 4x4 correlation matrices. For instance, 1 and 2 correspond to the combinations of (1,7,16,35) 



and (3,7,16,35) respectively and so on, as discussed in the text following Eq.(3.1). Masses are extracted 
according to the selection criteria described in Ref. [pi]]. The horizontal lines are drawn for the average 
mass over the four combinations (from 3"* to 6* as shown inside the shaded box). 



count is observed here as in Ref. [|n|] for a few of the interpolator basis smearing sets. However 
the ground and first excited states are robust against changes in the interpolator basis, providing 
evidence that an energy eigenstate has been isolated. It should be noted that the highest excited 
state (the third excited state) is influenced more by the level of smearing than the lower excited 
states. This is to be expected as this state must accommodate all remaining spectral strength and 
this is dependent on smearing. 

The P' combination in Fig.|l] provides heavier excited states as this basis begins with a low 
number of smearing sweeps (a sweep count of 1) and also contains another low smearing set of 
7 sweeps. The first and second excited states sit a little high in comparison with the other bases. 
This is also evident in Fig.^ for the lighter quark mass. Hence, extracting masses with this basis is 
not as reliable as other sets containing a great diversity of large numbers of smearing sweeps. The 
2"'^ combination also contains elements with a small smearing sweep count (3 and 7), hence this 
basis also provides heavier excited states and shows some systematic drift in the second excited 
state. However, this basis has reduced contamination from the excited states when compared with 
the first basis. 

We can observe at this point that including basis elements with 2 consecutive low smearing 
sweep counts (for instance, consecutive low numbers of smearing sweeps 1,7 of combination 
and 3,7 of 2"'* combination, respectively), provides a basis which does not span the space well. We 
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Figure 2: (Color online). As in Fig.|l|, but for the light pion mass of 249 MeV. 
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Figure 3: (Color online). The ground, the Roper and second excited states and the non-interacting P-wave 
N + n are illustrated. The black filled symbols are the experimental values of the ground and the Roper 
states. 
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also observe that the inclusion of basis elements with 2 consecutive high levels of smearing (for 
instance, a sweep count of 35, 48 as in the 1^^ combination) does not span the space well and gives 
rise to larger uncertainties. 

The 3'^'*, 4"^ and 5"^ combinations are well spread over the given range of smeaiing sweeps. 
They don't include successive lower smearing sweep counts. The 5"^ combination contains the 
basis element with a sweep count of 48 but has only slightly larger statistical errors than the 4"^ 
basis choice. All these bases provide diversity. It is observed that the 3'^'' through the 6* combina- 
tions provide consistent results for the first and second excited states. An analysis is performed to 
calculate the systematic errors associated with the choice of basis over the preferred four combina- 
tions (from 3'''^ to 6*) with O}, ■ 



j^j—^ {Mi — MY, where, Nb is the number of bases, in this 
case equal to 4. The statistical and systematic errors due to basis choices are added in quadrature. 
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Figure 4: (Color online). Compilation of current lattice results of the nucleon, , for the ground and first 
excited states reported in Refs. [| | |l [lOl 111 |23, |1 as described in the text. 



In Fig. ^ it is interesting to note that the observed lattice Roper state sits lower than the P-wave 
N indicative of attractive ttN interactions producing a resonance at physical quark masses. 

In Fig. ^ the ground state results are consistent for all the reported works. However, signifi- 
cant differences are readily observed in the results of the first positive parity excited (Roper) state. 
In Ref. [Q], Mathur etal. used a constrained curve fitting method. In Ref. [^, Sasaki etal. used the 
Maximum Entropy Method. In Ref. [19], Guadagnoli etal. used a modified correlator technique. 
Refs. [|[, ^ 21, 22, 23] all used the variational approach. Brommel etal. pO| ] considered 
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standard nucleon interpolators , Xi and and employ Jacobi smeared sources. They have per- 
formed simulations in a range of pion masses of 270-866 MeV and used 3x3 correlation matrices. 
They reported results that are too high to be interpreted as the Roper state. The early analysis of 
CSSM Collaboration by Melnitchouk etal. | [Io| ] used X\ and Xi interpolators and used a single 
Gaussian source smearing of 20 sweeps. They found no evidence for the Roper state as the energy 
states of their 2x2 correlation matrix analysis also proved too high to be interpreted as the Roper 
resonance. The results from Lasscock etal. from 3x3 correlation matrix analysis of standard 
interpolators also sits as high as those of Brommel and Melnitchouk, and cannot be interpreted 
as the Roper state. Burch etal. [22] considered the alternative approach of using Jacobi smeared 
sources of two different widths (narrow and wide) to increase the basis of operators and performed 
6x6 correlation matrix analysis for a pion mass down to 450 MeV. Though their lightest quark 
mass results sit slightly above our results labeled Mahbub etal. in Fig. |^ in the heavy quark- 
mass region the significant overlap between the results of Burch etal. and ours is apparent. Basak 
etal. 1 23] used non-local operators to form the basis of their correlation matrix and simulated at a 
pion mass of 490 MeV. Their results are also high. They reported that they did not find any Roper 
like positive-parity excitation. Using smeared-smeared correlators to construct a basis of correla- 
tion matrices, we identified a low-lying Roper state with significant curvature apparent at lighter 
quark masses [P]. 



5. Conclusions 



Through the use of a variety of smeared-smeared correlation functions in constructing corre- 
lation matrices, the first positive parity excited state of the nucleon , the Roper state, has been 
observed for the first time using the variational analysis |^. The current status of results for the 
Roper resonance from a number of groups are reviewed. Our lattice Roper state has a tendency to 
approach the physical Roper state showing a significant curvature as the chiral limit is approached 
[Q]. This work signifies the importance of using a diverse range of smeared-smeared correlation 
functions when constructing correlation matrices for the identification of the elusive Roper state. 
This robust approach should also be applied for larger dimensions of the correlation matrices not 
only built from the X\Xi correlators but also using X\X2 and in the negative parity channel. This 
will be the subject of future investigations. 
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